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ABSTRACT 

A new result of investigations of the hyperfine structure (hfs) anomaly in Be 
isotopes is presented. The hfs constant for n Be is obtained by using the core 
plus neutron type wave function: \2si > +\lds x 2 + ; | + >. A large hfs anomaly 
of 11 Be is found, which is mainly due to a large radius of the halo single particle 
state. 
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1 Introduction 



Recently, much interest has been paid to the magnetic hyperfme structure (hfs) 
for various nuclear isotopes H], |, H U, H El- For the experimental side, there is 
some progress in observing the transition of the hyperfme levels for atomic ground 
states. In fact, the accuracy of the measurement of the hfs splitting is improved 
a great deal, and even for lighter nuclei, there is a good chance of observing the 
hfs anomaly. This becomes possible due to the ion trap method which can isolate 
the atoms. This ion trap method [| can measure the hfs anomaly with the 
accuracy of order 10 -6 . 

The hyperfme structure has a sensitivity to the magnetization distribution 
in nucleus. It can, therefore, present a unique way to measure the neutron 
distribution in nucleus. In particular, there is a strong evidence that the n Li has 
a neutron 

halo[[| which may extend quite far over the typical nuclear radius of the 
neighboring light nuclei. Moreover, n Be has an anomalous spin-parity state for 
the ground state |l(J, and there may be some chance that it also has a large 
neutron radius. 

In this paper, we present a model calculation of the hfs anomaly for Be iso- 
topes. For 7 Be and 9 Be nuclei, we can use the Cohen-Kurath wave functions ||11|| . 
On the other hand, n Be has an anomalous spin-parity state, which is | + in- 
stead of | . A simple minded shell model wave function does not give a proper 
structure of the ground state. 

Here, we calculate the matrix element of the hfs operator for the ground state 
of n Be using the core plus neutron type wave function: a|2si > +{3\lds x 2 + : 
| + >. It turns out that the hfs anomaly e for n Be is quite large 

in magnitude compared to those of 7 Be and 9 Be. This is mainly related to 
the fact that the ground state of n Be has a large radius compared to those of 
ground states of other Be isotopes. 

The paper is organized as follows. In the next section, we briefly explain the 
theory of the magnetic hyperfme structure in electronic atoms. Section 3 treats 
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the core plus neutron type wave function. We evaluate the matrix elements of 
the hfs operator using this wave function. In section 4, numerical results of the 
isotope shifts of the hfs anomaly for Be isotopes are presented. Summary is given 
in section 5. 

2 Magnetic Hyperfine Structure 

An atomic electron which is bound by a nucleus feels the magnetic interaction 
in addition to the static Coulomb force. The magnetic interaction between the 
electron and the nucleus can be described as 

H ' = ~j JJv(r) • A(r)d 3 r (2.1) 
where the nuclear current jjv(r) can be written as 

Jiv(r) = ^-E gfV x Sl 5(r-R,)+£ (P^r - R.) + 5(r - R,)P,) (2.2) 

i i 

A(r) denotes the vector potential which is created by the atomic electron, and it 
can be written as 

A(r) = / j^dV (2.3) 
where j e (r) denotes the current density of an electron, and can be written as 

j e (r) = -ea<5(r-r e ) (2.4) 

where a denotes the Dirac matrix. 



2.1 The hfs anomaly 

The magnetic hyperfine splitting energy W can be written as 

W =<IJ : FF\H'\IJ : FF >= ^ [F(F + 1) - /(/ + 1) - J {J + 1)] aj (2.5) 

where I, J and F denote the spin of the nucleus, the spin of the atomic electron 
and the total spin of the atomic system, respectively, ai is called the hyperfine 
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structure (hfs) constant. Following ref. [E|, we can write the expression for the 
dj as 

a/ = 4 0) (l + e) (2.6) 
where of' is the hfs constant for the point charge, and can be written as 

4 0) = - r^V 1^ F^G^dr (2.7) 
1 J{J + 1) Jo 

where /i is the magnetic moment of the nucleus in units of the nuclear magneton 
fj,jsf, and F^ kJ ^ and G^ kJ ^ are the large and small 

components of the relativistic electron wave function for the kJ state, e is 
called hfs anomaly and can be written as 

0.626( feJ ) ^/i?A 2 0.386( fcJ ) * ( R,\ 2 (i)v( i )lrr 

(2.8) 

where is the single-particle operator of the magnetic moment, that is, fa = 
gi^Si + g\ £i, and is defined as 

E{ 1) =s i + V^[a* r < i,) ]j 1) . (2.9) 

i?o is a nuclear radius and can be given as i?o — with r = 1.2 fm. On the 
other hand, b^ kJ ^ is a constant which is calculated in terms of relativistic electron 
wave functions^], and is given as 

POO 

b {kJ) =0.23*^07(1 - 0.2 7 2 )/ / F {kJ) G {kJ) dr. (2.10) 

JO 

m e denotes the electron mass, k$ is a normalization constant, and 7 = Za. 
The isotope shift of the hfs anomalies of the two isotopes A 12 is defined as 

A 12 = ^-l (2.11) 
a h gi 

where g\ and g 2 are the total nuclear g-factors for isotope 1 and 2, respectively. 
Since the hfs anomaly e is quite small, A12 becomes 

A 12 ^ ei -e 2 . (2.12) 

The hfs anomaly e can be calculated once we know the nuclear wave function. 
Here, we employ simple-minded shell model wave functions with core polarization 



taken into account. We take the following approach which considers only the 
A£ = core polarization for the E^ operator. In this case, we can obtain the 
matrix element of the sj 1 ^ without introducing any free parameters as discussed 
in ref. 1112 . 



In this case, we can express the effect of the core polarization on the Ej operator 



we 



2.2 A£ = core polarization 

In this case, we can express the effect of 1 

in terms of the core polarization of the magnetic moment. Following ref. [12 
can write the expectation value of the E^ as 

< II\±^\II >= ±# K '|^+ 4( j!^' (ra) (M-^-^°-) (2.13) 

for I — £ ± ~. Here, g^ N ^ denotes the g-factor of the valence nucleon for the 
single particle state we are considering. fi sp is the single 

particle value of the magnetic moment. <5/i mes comes from the meson exchange 
current and can be given approximately as 

5^ mes « 0.1t 3 £. 

Therefore, we do not have any free parameters in the evaluation of the expectation 
value of the E^ . As the Be isotopes have orbits with small £, the exchange current 
effects are not important and we can safely neglect the term Sfj, mes . 



3 Core plus neutron model 



Recently, Suzuki et al.|T^] investigated the magnetic moment of the ground state 
of n Be. They describe the ground state of n Be as 

1+ 1+ 
| n Be I - j >= a| 10 Be(0 + ) x u2si : - > + / 3| 10 Be(2+) x v2ds : - > . (3.1) 



In this case, the magnetic moment of the Be ground state can be expressed as 



2 7 2 1 2 

\i = a fi v (28i) + V>v(ld B ) ~ qP ^( 2+ ) 

i 10 2 o 



(3.2) 



where fi v (2 Sl ), AVid s ) and fJ><2+) denote the magnetic moment of the neutron 2s i, 
Ids states and the 2 + (3.37 MeV) state of 10 Be, respectively. Here, we evaluate the 
hfs anomaly e for the two- component wave function, eq.(3.1). The hfs anomaly 
of eq.(2.8) can be written as 



0.62M ls ) 



OC /^(2si) < 



R 



R 



0.386 (ls) 
/i 



where 



2V ^ f R 



-J > 2s + T5 P < [jr) >u -3/3 W < ( 

(3.3) 



>-, • ./•/ Y,<hX U J.I > 



(3.4) 



and i runs over the nucleons in the state that has the total angular momentum j. 
MK2si) = l^v{id 6 ) — i#sj where = —3.826 is the spin g-factor for the neutron, 



and fi(2+) = 1.787JTJ]. The values of £„(2 Sl ), Ej,( lds ) and £( 2 +) are obtained as 



J ^(2si) 
2 



i/(ld s ) 
2 



2 2 

77<7$ n ), and £(2+) = —1.034, respectively. 



l 4 y s , — -(2+) : 

The H(2+) and the £(2+) are obtained by using the Cohen-Kurath wave func- 



tion, TBE(8-16)(TT|, 15]. Now, the recent measurement of the magnetic moment 
of n Be gives[|TJ[ 

n( n Be) = -1.682(3) fx N . (3.5) 
This indicates that the value of a 2 is close to 



a 



0.5 



(3.6) 



from the comparison of the observation with Fig. 2a of ref. |I4|. In this way, we 
obtain the hfs anomaly e for n Be. For the 7 Be and 9 Be isotopes, we can calculate 



the hfs anomaly using the Cohen-Kurath wave function [p4| . The hfs anomaly e 
for this case is written similarly as 



We also consider the case in which the effective g factor for the spin operator, 
g e Jf , is taken into account. As for the halo orbit, we assume that the second- 
order effects are rather small and result in little quenching of the spin operator. 
We, therefore, take g e J* = 1.0 for the z/2si/2- orbit. In this case, fi u (2s 1/2 ) = 
\gW = -1.913, /i Kld5/2) = igWgt" = ~lM3gf f , /i (2+) = 1.076 + 0.711 g e J f , 
and E„ (2 . 1/a) = \g^ = -1.913, E„ (1 * /a) = ^g e s ff = -2.460^, S (2+) = 
— l.OMg^f, in units of /in- The magnetic moment is, then, given by 

fj, = -1.913a 2 - (0.3585 + I.l298g e / f )f3 2 . (3.8) 

In comparison with the experimental value of the magnetic moment, —1.682(3) 
/ijv, we obtain a 2 = 61%. When the contributions from the meson exchnge cur- 
rents are taken into account, a 2 is close to 65%. We give results of the hfs anomaly 
e for 11 Be also for the case with the effective spin g factor. Numerical evaluations 
are carried out in section 4. 



4 Numerical Results 

Now, we evaluate numerically the hfs anomaly e for the Be isotopes. 

4.1 Core plus neutron type wave function 

In order to evaluate eq.(3.3), we need to know the magnetic moment \i for the Be 
isotopes. The magnetic moments of 9 Be and n Be are observed, but the magnetic 
moment of 7 Be is not yet determined experimentally. We evaluate the magnetic 
moment of 7 Be empirically. 

For 7 Be, we can make use of the magnetic moment of 7 Li since they are isospin 
doublet states. In this case, we predict the magnetic moment for 7 Be 

/i Pred ( 7 Be) = -1.377. (4.1) 
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In Table 1, we list the values of the quantities which are necessary to calculate 
the hfs anomaly e. The r.m.s. radii are those for wave functions solved in a 
Woods-Saxon potential with Rq = 1.2 A 1 / 3 and a=0.60 fm. For the 2si/2 and 
M5/2 neutron orbits in n Be, the wave functions are obtained to reproduce the 
separation energies, 0.50 MeV and 3.87 MeV, respectively. The anomaly e for 
n Be is obtained from eq. (3.3) to be 

e( n Be) = -0.12015a 2 - 0.02326/3 2 (%). (4.2) 

It gives e = —0.0717% for a 2 = 0.50 and [3 2 = 0.50. Magnitude of e gets as large 
as -0.091% (-0.101%) as a 2 becomes 0.70 (0.80). The hfs anomalies for 7 Be and 



9 Be are obtained from eq.(3.7) by using the Cohen-Kurath wave functions |TT|, |i~5|j . 
In Table 2a, we present the calculated values of hfs anomay e for the Be 

isotopes. Now, it turns out that the hfs anomaly for 11 Be has a very large value 
compared to other isotopes. This is mainly connected with the fact that the 
n Be has a large neutron radius. The r.m.s. radius of the halo 2si orbit 



becomes as large as 6.4 fm in deformed Woods-Saxon models^]. When we use 
a value of 6.37 fm for the r.m.s. radius of the halo 2s 1 orbit, which is obtained 



in the deformed Woods-Saxon model^, [TJ], the hfs anomaly is given by 

e = -0.12827a 2 - 0.02326/3 2 (%). (4.3) 

Eq. (4.3) leads to e = —0.0758 % for a 2 = 0.5, which is close to the value obtained 
from eq. (4.2). The anomaly e becomes —0.076 ~ —0.097 % for a 2 = 0.5 ~ 0.7. 
We finally give numerical results for the case with the effective spin g factor. We 
take g e J* = 1.0 for the halo z^2si/2 orbit and g e J? = 0.85[pl|] for other normal 
orbits. The hfs anomaly is given by 

e = -0.12015a 2 - 0.020 14/3 2 (%). (4.4) 

Here, the 2si/2 orbit obtained in the spherical Woods-Saxon potential 

is used. The anomaly becomes e = —0.0852% for a 2 = 0.65, whose magnitude 
is larger than the value obtained from eq. (4.2) about by 0.01 %. The anomaly 
e becomes -0.070 0.100% for a 2 = 0.5- 0.7. 
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4.2 Single particle model with core polarization 



Next, we calculate the hfs anomaly e for Be isotopes using the single particle 
model with core polarization. Since the radius of the single particle state can be 
different from that of 

the core polarization state, we modify the expression eq. (2.8) in the same 
way as eq. (3.3) 

R\ 2 „ fR^ 2 



0.626( ls ) 



0.386 (ls) 



<yN )Ti 



up 



< 



R^ 2 

Rq 



>s P +^5hcp < 



where S sp denotes the expectation value of < II \ J2t=i ^ 
single particle state and can be written as 

v __3(/ + |) „ . 1 



R_ 

Rq 

(1) 



>CP 



(4.5) 



for / 



+ 



II > with the 



(4.6) 



bp 4(7 + 1) — ~ ' ' 2" 

S/j,cp is the magnetic moment which arizes from the core polarization f20|. < 
(-^=r) 2 >cp denotes the expectation values with the states involved in the core 
polarization. In Be isotopes, they are lp states. In Table 2b, we list the calculated 
values of the hfs anomaly e for the Be isotopes. In the same way as the core 
plus neutron type calculation, the hfs anomaly for 11 Be has a very large value 
compared to other isotopes. This is due to the fact that the n Be has a large 
neutron radius since it has an anomalous spin-parity state. 



5 Conclusions 

We have presented the numerical calculations of the magnetic hfs anomaly for the 
Be isotopes. First, we employ the wave function which has a component coupled 
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to the 2 + — core excitation. This gives a large hfs anomaly for 11 Be since the 
neutron outside the shell is assumed to be 2s 1 or Ids orbits. On the other hand, 

2 2 

we find that the 7 Be and 9 Be isotopes have a small hfs anomaly. We also evaluate 
the hfs anomaly using a single particle shell model with core polarization. We 
also predict a very large hfs anomaly for n Be. 

It would be extremely interesting to learn whether the very large hfs anomaly 
of n Be can be realized in nature or not. Since this is related to the radius of the 
neutron halo nucleus, it may well help understand the structure of the neutron- 
rich nuclei. We hope that experimental observations will clarify this point in near 
future. 
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Table 1 

Ro, b*- 18 -*, r.m.s. radius and magnetic moment 





7 Be 


9 Be 


n Be 


i? [fm] 


2.296 


2.496 


2.669 


[%} 


0.0170 


0.0185 


0.0198 


< r2 >(ip) [H 


2.553 


2.569 


2.588 


< r 2 >f 2s) [fm] 






6.165 


< r2 >h) H 






3.551 


f^exp 


n.m.] 


-1.377* 


-1.177 


-1.682 




n.m.] 


-1.913 


-1.913 


-1.913 




n.m.] 


0.536* 


0.736 


0.231 



The values of R , b^ ls \ r.m.s. radii for the lp and 2s states and the 
magnetic moments for Be isotopes are shown. * indicates that the value 
of the magnetic moment is empirically extracted from those of other 
nuclear isotopes. 



Table 2a 

hfs anomaly e and isotope shift A i2 





7 Be 


9 Be 


n Be 


e [% 




-0.0245 


-0.0249 


-0.0717 


&7,A [ 


%} 





0.0004 


0.0472 



The calculated values of the hfs anomaly e and the isotope shift A 12 for 
the Be isotopes obtained from eqs. (3.3) and (3.7) are shown, a 2 = j3 2 
= 0.50 are used for n Be. 



Table 2b 

hfs anomaly e and isotope shift A 12 
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7 Be 


9 Be 


n Be 


e [% 




-0.024 


-0.023 


-0.118 


&7,A [ 


%] 





-0.001 


0.094 



The calculated values of the hfs anomaly e and the isotope shift A 12 
the Be isotopes obtained from eq. (4.5) are shown. 
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